Dynamic Correlation in Wave Propagation in Random Media 
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We report time- resolved measurements of the statistics of pulsed transmission through quasi-one- 
dimensional dielectric media with static disorder. The normalized intensity correlation function 
with displacement and polarization rotation for an incident pulse of linewidth a at delay time t 
is a function only of the field correlation function, which is identical to that found for steady- 
state excitation, and of Ha{t), the residual degree of intensity correlation at points at which the 
field correlation function vanishes. The dynamic probability distribution of normalized intensity 
depends only upon Ko-(t). Steady-state statistics are recovered in the limit ct^O, in which «;o-=o is 
the steady-state degree of correlation. 
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The statistics of steady-state classical wave propaga- 
tion Q, 1^ 1^ and electronic conductance 4] in disordered 
systems reflect the superposition of partial waves follow- 
ing trajectories with a wide distribution of path lengths, 
which is proportional to the particle time-of-flight dis- 
tribution. Since the lengths of trajectories increase with 
time, paths cross upon themselves more frequently and 
the impact of weak localization could be expected to 
build in time dllll^, 9, 10, 11, H E [11 [3. This 
was demonstrated recently in the observation of an in- 
crease of the coherent backscattering enhancement from 
2 to 3 in time-resolved acoustic measurements in a three- 
dimensional elastic body 9] , and in a time-decaying leak- 
age rate of microwave radiation from a quasi-lD ran- 
dom dielectric sample 12]. To achieve a systematic un- 
derstanding of weak localization in the time domain, it 
is essential to examine the statistics of propagation in 
addition to ensemble- averaged transport. This can be 
accomplished by parsing transmission according to the 
delay from an exciting pulse and studying the correla- 
tion and probability distribution of intensity as a func- 
tion of delay time, and finally by relating these dynamics 
to steady-state statistics obtained under monochromatic 
excitation. A dynamical perspective on the impact of 
gain or loss on steady-state statistics might be of partic- 
ular interest since the distribution of trajectories within 
the sample at a given delay is not altered by the presence 
of inelastic processes, while the average temporal profile 
changes. 

Nonlocal intensity correlation iH ll^ ll^ fl^ leads to 

f it transmission fluctuations 3l lid Eol 
liM ii^l and lies at the heart of mesoscopic physics 
. The intensity correlation function associated with 
bulk scattering versus displacement or polarization rota- 
tion of a detector may be expressed as C=F+k{F + 1) 
[II III Ei, where F^Fe]'^ is the square of the field 
correlation function. The degree of correlation, k, is the 
value of C when F vanishes, which occurs, for exam- 
ple, for a polarization shift of 90° or for displacements 



much greater than a wavelength. In the absence of in- 
elastic processes, the probability distributions of both 
intensity and total transmission normalized to their re- 
spective ensemble averages, P{sab) and P{sa), where a 
and b are modes of the incident and transmitted waves, 
are obtained from diagrammatic 20.] and random matrix 
|2l[ theories in terms of the dimensionless conductance, 
g !25|, in the diffusive limit, (?2>1. For both quantum 
and classical waves, g is equal to the sum of transmission 
coefficients over all input and output modes (2dj . Sur- 
prisingly, theoretical expressions for P{sab) and P (s„) 
\2(i \2l\ closely match the measured distributions [23] 
even in the presence of absorption, and even at the local- 
ization threshold, reached when g^l, when g is replaced 
by 2/3var{sa)- 

In this Letter, we report microwave measurements 
of the time-resolved field transmitted through random 
quasi-lD dielectric samples. Remarkably, the field corre- 
lation functions with displacement and polarization ro- 
tation at any time are identical to those found in steady 
state. We also find that the corresponding cumulant cor- 
relation functions of the normalized intensity have the 
same dependence upon the field correlation functions as 
in the frequency domain, Ca{t)=F+Kc{t){F + 1), with 
a parameter, Ka{t), expressing the degree of correlation 
at time t following excitation by a Gaussian pulse of 
width (T. We find further that the probability distribu- 
tion of normalized intensity, P{sab{t)), has the same form 
as for the steady-state transmitted intensity distribution 
[23,123, but with 2/3K<,(t) substituted for g. We find 
that even in diffusive samples, Ka{t) reaches values ex- 
ceeding the steady-state value at the Anderson localiza- 
tion threshold [23 of k~2/3 [2^. Steady-state statistics 
are found to be a limiting case of dynamic statistics, in 
which the incident pulse linewidth vanishes, a—^0, with 
Ka-=o=K. These results show that Ka{t) is the essential 
function describing the statistics of wave propagation. 

Spectral measurements of the field transmission coeffi- 
cient of microwave radiation as a function of rotation of 
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linear polarization and displacement are made with use of 
a vector network analyzer in ensembles of random dielec- 
tric samples in which the wave is diffusive. The samples 
are contained in a copper tube with open ends of length 
L greatly exceeding its 7.3-cm diameter. New sample 
realizations are produced by briefly rotating the sample 
tube about its axis after each spectrum is taken. The 
response to a pulse with a Gaussian temporal envelope 
at carrier frequency is obtained by Fourier transform- 
ing the product of the field transmission spectrum and a 
Gaussian spectral function of width a. 

Polarization-selective measurements of the transmitted 
microwave field are made with use of a conical horn detec- 
tor The samples are composed of 0.95-cm-diameter 
alumina spheres with refractive index 3.14 embedded in 
Styrofoam spheres of diameter 1.9 cm and refractive in- 
dex 1.04 at an alumina volume fraction of 0.068. Mea- 
surements are made in an ensemble of 10,000 alumina 
configurations of L=61 cm (Sample A) over the frequency 
range 14.7-15.7 GHz in steps of 0.6 MHz for a single ori- 
entation of the horn detector. Measurements are also 
carried out in an ensemble of 12,000 alumina configura- 
tions of L=90 cm (Sample B) over the frequency range 
16.95-17.05 GHz in 1-MHz steps, for 7 orientations of the 
horn detector rotated in steps of 15° over a 90° range. 
Steady-state measurements of intensity correlation give 
K=0.09 in Sample A [13 and /t=0.29 in Sample B [ig |. 

Examples of the response to excitation pulses with 
three different values of cr in a single random realiza- 
tion of Sample A are shown in Fig. la versus time de- 
lay from the center of the incident pulses in units of 
the diffusion time, t o={L+ zi^Y / tt"^ D , where D is the 
diffusion coefficient and zq is the boundary extrapola- 
tion length. The width of temporal intensity fluctua- 
tions is seen to be approximately equal to that of the 
incident pulse. This can be expressed quantitatively 
via the correlation function with time shift At of the 
transmitted field normalized to the square root of the 
ensemble average of the time-varying intensity at time 
t, Fe„ = {E* {t)E^ (t+At)) /V {la (t)) iC (t+At)) . We find 
that Fe„ is independent of t. In the case of a Gaussian 
incident pulse of bandwidth a, \Fe^ \ is a Gaussian with 
width {V2TTa)-^ (Fig. lb). We also find that, though 
la-it) in any given realization depends strongly upon a, 
(la-it)) for t > t]j depends only weakly upon cr, once 
cr > l/n'^tD- 

The correlation function of the normalized transmit- 
ted field with shift in polarization angle A9, Fe(A0), is 
shown in Fig. 2a for the two delay times indicated by 
the vertical lines in the inset for the incident pulse of 
bandwidth ct=5 MHz and for monochromatic excitation. 
All functions are well described by the steady-state re- 
suh, Fe(A6')=cos(A6') UlEil. The corresponding in- 
tensity correlation functions are shown in Fig. 2b. These 
have the form, C^(A6l, t)=F(Ae')+K^(i)[F(A6') + l], with 

Kait)=Cai90,t). 
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FIG. 1: (a) Transmitted intensity through a random realiza- 
tion of Sample A (solid curves) following incident Gaussian 
pulses (dashed curves) with cr=7, 15, and 30 MHz. The in- 
cident pulses are centered at t=0 and shown with the same 
height, (b) Absolute value of the field correlation function of 
Ea (t) /-\/ (la-it)) with time shift for the values of a in (a) . 



In order to study the dynamics of spatial correlation, 
measurements are taken at 50 points separated by 1.06 
mm along a line centered on the tube axis on the output 
surface of a random sample of Polystyrene spheres Il8| . 
A 3-mm antenna is aligned perpendicular to the line of 
displacement. The Polystyrene spheres of refractive in- 
dex 1.59 are packed at a volume filling fraction of 0.52 in 
a copper tube with i=100 cm. Measurements are made 
in 1380 sample realizations over the frequency range 17.2- 
17.8 GHz in steps of 0.625 MHz. In this sample, k=0.06. 

The correlation function of the normalized trans- 
mitted field with displacement Ar, ^^(Ar), is shown 
in Fig. 3a for the two delay times indicated in 
the inset following pulsed excitation with cr=20 MHz 
and for monochromatic excitation. The overlap 
of these curves indicates that the dynamic FsiAr) 
is independent of t and identical with the steady- 
state field correlation function. The correspond- 
ing intensity correlation functions, shown in Fig. 3b, 
have the form, C„(Ar, i)=F(Ar)+K<,(t)[F(Ar)-f 1], 
with Kait) determined from the residual correla- 
tion at displacements Ar > 3.5 cm, except for 
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FIG. 2: (a) Real part of the field correlation function, 
Fe{A6), and (b) intensity correlation function, Co-(A6), 
with polarization rotation of the transmitted wave through 
Sample B at the two delay times following pulsed exci- 
tation with a— 5 MHz and for monochromatic excitation 
(CW). The dashed curves are (a) Fe(A6I)=cos(A6') and (b) 
C^{Ae,t)=F{Ae)+Ka{t)[F{Ae)+l] with the values of K^it) 
indicated. The solid curve is F{A9). The logarithm of the 
average pulsed transmission through Sample B for a— 5 MHz, 
normalized by the average steady-state transmitted intensity, 
is shown in the inset. The delay times at which correlation is 
presented in the figure are indicated by vertical dashed lines. 



FIG. 3: (a) Real part of the field correlation function FsiAr), 
and (b) intensity correlation function, Ca{Ar), with displace- 
ment on the output surface of a random Polystyrene sam- 
ple at the two delay times following pulsed excitation with 
a— 20 MHz and for monochromatic excitation (CW). The 
dashed curves in (b) are C:y{Ar,t)^F{Ar)+n^{t)[F{Ar) + l], 
with the values of Kcr{t) indicated. The solid curve is F{Ar). 
The logarithm of the average pulsed transmission through the 
sample for a— 20 MHz, normalized by the average steady-state 
transmitted intensity, is shown in the inset. The delay times 
at which correlation is presented in the figure are indicated 
by vertical dashed lines. 



t/irf=4.8 when Ka{t) is determined from the relation, 
var(safc(t))a=a(0, t)=l+2K^{t). 

The time-resolved probability distributions of normal- 
ized transmitted intensity, P{sab{t)), are shown in Fig. 4 
for the three samples considered, for various values of de- 
lay time and pulse bandwidth. The steady-state inten- 
sity distributions for Sample B and for the Polystyrene 
sample are also shown. Apart from the uppermost curve, 
each of the curves is displaced by a multiple of one decade 
for clarity of presentation. The dashed curves are ob- 
tained from the expressions for the steady-state intensity 
distribution P{sab) 20, 21] in the limit, 5^1, and in the 
absence of absorption, but with 2/3Kcr substituted for g, 

P{Sab) = / -P{Sa)exp{~Sab/Sa) , (1) 



with 

/ioo 1 
— exp[us,-<i>(u)], (2) 
-ioo 27ri 

where 

= (2/3k^) \n^(y/l + ?,vK„/2 + . (3) 

The values of Ko- for Sample B and for the Polystyrene 
sample are obtained from the measured intensity cor- 
relation functions in Figs. 2 and 3. The values of 
K„ for Sample A are determined from the relation, 
Ca{Q,t)=l+2Hicr{t). Excellent agreement with measured 
results is obtained in all cases. These encompass both 
steady state and dynamic propagation, in the presence 
of strong and weak absorption, in the weak as well as the 
strong correlation regimes, in which exceeds its value 
at the Anderson localization threshold of 2/3. 
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FIG. 4: Dynamic and steady-state (CW) probability distribu- 
tions of normalized transmitted intensity (solid curves) with 
the values of t/to, cr, and Ko-, and sample type indicated. 
The dashed curves are from Eqs. (1-3). In the case of distri- 
bution (1) of long delay time, noise accounts for nearly half 
the measured field. The curve through the data is obtained 
by transforming accordingly Eqs. (1-3) to include the effect 
of noise. In the cases other than (1), noise is disregarded. 



In conclusion, the time-resolved intensity correlation 
function has a universal form in terms of K^{t) and the 
square of the time-independent field correlation function, 
F. At the same time, the probability distribution of in- 
tensity is determined exclusively by Ka-(t). In the limit, 
CT^O, Ka~^K, yielding steady-state statistics. Thus, the 
time- varying degree of correlation, K^{t), is the control- 
ling function of mesoscopic statistics. It is, therefore, of 
prime importance to explore the possibility of a universal 
formulation of the time variation of the degree of corre- 
lation and its relationship to spatial localization. 
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